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ABSTRACT 

Pulsars are remarkable objects that emit across the entire electromagnetic 
spectrum, providing a powerful probe of the interstellar medium. In this study, 
we investigate the relation between dispersion measure (DM) and X-ray absorp- 
tion column density A^h using 68 radio pulsars detected at X-ray energies with 
the Chandra X-ray Observatory or XMM-Newton. We find a best-fit empirical 
hnear relation of A^h (10^°cm~^) = 0.30iQ;o9 DM (pccm~^), which corresponds 
to an average ionization of 10^3%, confirming the ratio of one free electron per 
ten neutral hydrogen atoms commonly assumed in the literature. We also com- 
pare different iVn estimates and note that some TVh values obtained from X-ray 
observations are higher than the total Galactic Hi column density along the same 
hue of sight, while the optical extinction generally gives the best A^h predictions. 

Subject headings: dust, extinction — ISM: general — pulsars: general — X-rays: 
ISM 
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1. INTRODUCTION 

The broadband emission of pulsars from radio frequencies to 7-rays can be used to 
probe the physical conditions of the interstellar medium (ISM). Specifically, their radio 
pulsations allow accurate measurements of the free electron column density and their X-ray 
extinction traces the interstellar gas along the line of sight. Radio waves travelling in the 
ISM are dispersed by free electrons such that signals at lower frequencies propagate at a 
lower speed and hence arrive on Earth later than those at higher frequencies. The time 
delay (At) between two observing frequencies (z^i, 1^2) depends on the dispersion measure 
(DM), which is the integrated free electron number density Ug from Earth to the source at 
distance d: 

where nie and e are electron mass and charge, respectively, and c is the speed of light. 
Most free electrons in our Galaxy are found in the hot phase of the ISM, including Hll 
regions ionized by UV radiation from hot O or B type stars and the shock-heated interior 
of supernova remnants (SNRs). These sources can contribute significant DM up to a few 
hundred parsecs per cubic centimeter. At X-ray energies, photons are absorbed mostly by 
heavy elements in the interstellar gas due to the photoelectric effect. This has a strong 
energy dependence and is most prominent in the soft X-ray band. As a result, it modifies 
the observed low-energy portion of the X-ray spectrum and has to be accounted for in 
spectral modeling. The amount of extinction, which is expressed in terms of the equivalent 
atomic hydrogen column density A^h, is se nsitive to gas and molecular clouds, which traces 



the warm and cold phases of the ISM (see lWilms et al.ll2000l ). 



One natural question to ask is whether there is any correlation between DM and Ah 
in our Galaxy. Such a correlation can reflect the physical connection between different 
phases of the ISM. Also, it can provide a useful tool to estimate one quantity from the 
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other, help plan new observations and determine X-ray luminosity upper limits in cases 
of non-detection. In the literature, an average ionization fraction of 10% in the ISM, i.e. 
one free electron per 10 equ ivalent hydrogen atoms, has been comnaonly a ssumed in orde r 



to infer A^h from D M (e.g. 



Camilo et al. 



Seward fc Wang 



1988 



Kargaltsev et al. 



2007 



Gil et al. 



2008 



20121 ). but the justification for this choice has been unclear. X-ray-emitting 



radio pulsars offer a powerful diagnostic tool for a quantitative study of the correlation. 
Because they are model-independent and relatively straightforward to measure from radio 
timing, DM values are well determined, typically to better than a fractional uncertainty of 
10^'^. However, what has made the determination of any DM-A^h correlation difficult in the 
past is the lack of high-quality X-ray data for A^h measurements. In particular, previous 
generations of X-ray telescopes had poor angular resolution that precluded discerning the 
pulsar emission from that of the surrounding SNRs and pulsar wind nebulae (PWNe). 
Thanks to new X-ray missions such as the Chandra X-ray Observatory and XMM-Newton, 
precise measurements of A^h have been obtained for many pulsars in recent years, allowing 
a statistical study of A^h values for the first time. 

In this paper, we compile a list of DM and A^h values for 68 X-ray-emitting radio pulsars 
using the latest Chandra and XMM-Newton measurements reported in the literature. 
We found a clear correlation between these two column densities and obtained a best-fit 
empirical relation of 10^3% ionization. In Section 2, we describe our sample selection 
criteria. The statistical analysis and results are presented in Section 3, and we discuss the 
implications of our results in Section 4. 



SAMPLE SELECTION 



We started with a list o f X-ray detect e d rad io pulsars from iPossenti et al.l (120021 ). 



Becker fc Aschenbachl (|2002|), 



Pavlov et al. 



(120071 ) 



Kargaltsev fc Pavlovl (120081 ). and 
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Kargaltsev fc Pavlovl (120101 ). then expanded the sample through careful literature searches 
for updated observational results and recent discov e ries. The latter in clude three magnetars 



that show radio emission (jCamilo et al 



20061 . 



20071 : 



Levin et al 



2010l ) and over a dozen new 



pulsars identified in 7-rays with the Fermi Gamma-ray Space Telescope and subsequently 



Marelli et al 



2OI1I ). Finally to 



detected in follow-up radio and X-ray observations (see 
complete the list, we went through the Chandra and XMM-Newton data archive to search 
for pulsar observations, and looked up relevant publications based on these data. 



20051). 



The pulsar DMs are adopted from the ATNF Pulsar Catalog^ (iManchester et al 
They are all very well measured with negligible uncertainties compared to those for A^h- On 
the other hand, it is much more difficult to determine A^h, because this requires a strong 
X-ray source and good knowledge of the intrinsic emission spectrum. The X-ray emission 
of pulsars is not fully understood; commonly used models include a blackbody (BB) and 
a neutron-star hydrogen atmosphere (NSA) for the thermal emission, and a power law 
(PL) for the non-thermal emission. More complicated models consisting of thermal and 
non-thermal components are sometimes used. To minimize any bias, we selected the A^h 
values for our sample according to the following criteria: 



1. We restricted our choices to those in the latest studies using the Chandra and 
XMM-Newton observations, since the good angular resolution and sensitivity of 
these telescopes offer high-quality spectra with minimal background contamination. 
Any joint fits with other X-ray telescopes are not considered, in order to avoid 
cross-calibration uncertainties. 

2. We adopted only A^h values from actual X-ray spectral fits in which the A^h is allowed 
to vary freely, and ignored any A^h inferred from DM, optical extinction (Ay), or total 

^|http : / /www . atnf . csiro . au/research/pulsar/psrcat/ 
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Galactic Hi column density. 

3. A^H from the best-fit spectral model is always preferred, unless there are physical 
arguments favoring another model. If different emission models give the same 
goodness-of-fit and the authors do not indicate a clear preference, we choose the 
simpler one. For example, we prefer a BB model over an NSA model, since the latter 
requires more assumptions, including the atmosphere composition, surface magnetic 
field and gravity. 

4. For pulsars associated with bright PWNe, the nebular A^h values are adopted if they 
are better constrained than those of the pulsars, because the simple PL spectra of 
PWNe can reduce systematic uncertainties in spectral modeling. A^h from SNRs 
are used in a few cases when the pulsars and PWNe are too faint for useful A^h 
measurements. 

Our final sample contains 68 pulsars. One of them (PSR B0540— 69) is extragalactic 
and only two (PSRs J1740— 5340 and B1821— 24) are in globular clusters; cluster pulsars 
are generally too faint for precise A^h measurements. The pulsar DM and A^h values are 
listed in Table [1] and plotted in Figure [1] The reported statistical uncertainties and upper 
limits for A^h are at 90% confidence level, i.e. 1.6a. We list in the Table the X-ray spectral 
models used to obtain A^h- The choice of spectral model is clear in all cases except PSRs 
J1622— 4950 and B1757— 24, for which both thermal and non-thermal fits are acceptable. 
Nonetheless, A^h from different fits only varies by a factor of 2 for J1622— 4950 and does 
not change for B1757— 24. Therefore, we conclude that systematic bias induced by spectral 
models is minimal. 

Table [1] also shows the pulsar Galactic coordinates [l, b) and distances, and this 
information was used to calculate the vertical height {z) from the Galactic Plane. The 
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coordinates are taken from the ATNF Pulsar Catalog and distance estimates are obtained 
from parallax measurements, Hi absorption measurements of the pulsars or the associat ed 



SNRs, or DM using the NE2001 Galactic electron density model (ICordes Sz Lazio 



20021 ) ■ 



If available, parallax distances are always p referred since they ar e the most accurate. All 



parallax and Hi distances are adopted from 



Verbiest et al. 



(I2OI2I ) and references therein. 



and have been corrected for the Lutz-Kelker bias, e xcept for PSR J1023+0038, which has a 



recent parallax measurement by 



Deller et al 



(|2012| ). For DM distances, we did not attempt 
to derive the uncertainties, b ut note that the fractional uncertainties could be 25% or larger 
(see e.g., ICamilo et al.ll2009l ). Finally, there are exceptional cases in which previous studies 



argue for different distances than the DM-estimated ones. They are noted in the Table. 
The pulsar A^h and DM are plotted against distance in Figures [2] and [3l respectively. 



3. ANALYSIS & RESULTS 

Figure [1] shows a positive correlation between the pulsar DM and Ah values, with 
deviations ranging from a factor of a few to an order of magnitude. There are some obvious 
outliers, including the Vela pulsar (PSR B0833— 45), the double pulsar (PSR J0737— 3039), 
and PSR J1747— 2809 in the Galactic Center direction. To quantify the DM- Ah correlation, 
we ignored pulsars with Ah upper limits and obtained a Pearson's correlation coefficient 
of 0.72. This is significant since the one-tailed probability of such a correlation arising 
by chance from unrelated variables is only 4 x 10~^. More useful is an empirical relation 
between these two observables. We performed a linear fit to the data by minimizing 
the value. Ah measurements with fractional uncertainties larger than 80% or upper 
limits only (gray points in Figure [T]) are excluded in the fit. We also ignored the Vela 
pulsar, which is located in the Gum Nebula inside the hot and low-density Local Bubble, 
and PSR B0540— 69, which is in the Large Magellanic Cloud (LMC), because they seem 
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unlikely to follow the DM-A^h correlation as would other Galactic sources. Only statistical 
uncertainties in A^h are considered in the x^"fit since uncertainties in DM are negligible. 
Also, we did not attempt to model the systematic uncertainties, but we note that the ones 
introduced by different photoelectric absorption models and elemental abundances, or by 



cross-calibratioii betw een telescopes are only at a few percent level (see 



Wilms et al. 



2000 



Tsujimoto et al. 



20111 ). relatively small compared to the statistical uncertainties. Assuming 



A^H and DM are directly proportional, the best fit gives 

A^H (10^° cm-2) = 0.30+°:^^ DM (pc cm'^^) , (2) 
corresponding to an average ionization of 101^3%. "^^^ 90% confidence int erval is quoted here. 



which is obtained from 10000 simulations via bootstrapping resampling (lEfron &: Tibshirani 
19931 ). The result is plotted in Figure [H We also tried fitting a more general linear relation 
by fitting the y-intercept as well, but found that the latter is consistent with zero at 90% 
confidence. If we ignore the measurement uncertainties in A^h and perform a least squares 
fit, we obtain A^h (lO^'^cm"^) = 0.83 DM (pccm~^), giving a lower average ionization of 4%. 

To check if the DM-A^h relation could depend on the source location in the Galaxy, we 
divided the sample into groups according to their vertical height from the plane and their 
Galactic longitudes. The results are shown in Figures [T](b) and[T](c), respectively. In the 
high-DM regime, sources toward the Galactic Center direction, e.g., PSRs J1747— 2958 and 
J1747— 2809, show a hint of a larger A^n-to-DM ratio. However, the systematic variation is 
less clear at lower DM and our limited sample precludes a detailed analysis. In Figure [2] 
we plotted A^h against distance. This indicates a general correlation, albeit with a large 
scatter. There is also a hint that for sources at a similar distance, A^h is systematically 
larger near the Galactic Plane (Figure [SJ^b)), however, the dependency on Galactic longitude 
is less clear (Figure [2](c)). The DM variation with distance is presented in Figure [31 While 
this may seem to exhibit a good correlation at large distances, we note that sources with 
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DM-derived distances provide no new information, only the NE2001 model prediction. 
In addition, there is a very large range of DMs for nearby pulsars around 300 pc, from 
2.4 ± 0.2pccm^'^ for PSR J0108— 1431 to 68 ± 1.6pccm~'^ for the Vela pulsar, spanning 
nearly a factor of 30. Similar to A^h, Figure |3]^b) also indicates a higher DM toward the 
Galactic Plane. 



4. DISCUSSION 



We have investigated the DM-A^h connection for 68 radio pulsars detected with 
Chandra or XMM-Newton. We found a good correlation between these two column 
densities, suggesting that free electrons in the Galaxy generally trace the interstellar gas. 
That said, some A^h values in Figure [1] show significant deviation from the best-fit line, by 
a factor of a few up to an order of magnitude. This could be attributed to inhomogeneity 
of the ISM, possibly due to molecular clouds, supernova remnants, or Hll regions in the 
line of sight. Such an effect is more prominent for nearby sources, since the distribution 
of free electrons and interstellar gas is highly anisotropic around the Local Bubble (see 



Taylor fc Cordes 



1993 



Lallement et 



al. 



contribution from the Gum Nebula ^Taylor fc Cordes 



20031). In particula r, there is significant DM 



19931 ). resulting in a wide range of 



DMs for pulsars within ~300pc (e.g., the Vela pulsar and PSR J0737— 3039; see Figure [3]). 
At large distances, local fluctuations are expected to average out and the scatter of A^h 
and DM with respect to distance likely arises from Galactic structure, such as the disk. 



Cox 



20051). We 



spiral arms, and different scale heights of various ISM components (see 
have attempted to identify any systematic trends in DM and A^h with respect to source 
location. While Figures I21^b) and EJ^b) hint at higher A^h and DM toward the Galactic 
Plane, more sources are needed for a quantitative comparison with the detailed Galactic 
structure. Beyond our Galaxy, we note that while PSR B0540— 69 in the LMC was not used 
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in the fit, its DM-to-A^H ratio lies close to the best-fit line in Figure [H This is somewhat 
surprising because of th e different interstellar abundances in the LMC than in our Galaxy 



(IRussell fc Dopita 



19921 ). We argue that this could merely be a coincidence rather than 



the general 



(jPark et al 



case. Indeed, the LMC contributes 9 0% of the iVn toward P SR B0540-69 



2OIOI ) but only two thirds of the DM (iManchester et al 



20061) 



The DM-A^H correlation can be used to estimate one quantity from the other, offering 



a useful tool for pulsar observations. For instance, radio pulsations have 
the magnetar 4U 0142+61 with a DM of 27 ± 5 pc crn "^ (iMalofeev et al 



3een claimed from 



A^H value of 9.6 ± 0.2 x lO^^ cm~2 f lGohler et al 



20101 ). Given its 



2OO5I ) . the claimed DM seems somewhat 



small when compared to other sources of similar A^h in Figure [T] For X-ray observations, 
there are many cases requiring a priori knowledge of A^h, including flux estimates when 
planning for new observations, measuring the intrinsic spectra of faint sources, and deriving 
luminosity limits for non-detection. In many previous studies, A^h i s inferred from the DM 



by assuming one fr ee electron per ten neutral hydrogen atoms (e.g. 



Camilo et al. 



Kargaltsev et al. 



20071 : 



20121 ). Our result directly conflrms that this is a reasonable approximation. 



but as a caveat, the scatter in A^h is typically a factor of a few up to an order of magnitude. 



I n addition to DM, t he total Galactic Hi column density from 21-cm radio surveys 



(e.g., 
(e.g., 



Kalberla et al. 



Olausen et al. 



20051) and Ay have also been used as proxies for the X-ray absorption 



20131 ). These A^h estimates are plotted in Figure HI It is clear that 



some X-ray-inferred A^h values exceed the total Hi column density of the Galaxy. As 
shown in the Figure, the latter saturates at ~ lO^^cm"^, resulting in gross underestimates 
for high-DM (> lOOpccm"^) or distant (> 3kpc) pulsars. It has been reported that at 
high Galactic column densities > lO^^cm"^, which occur at low Galactic latitudes, the 



X-ray absorption columns are generally larger than the 



(lArabadjis fc Bregman 



1999 



il co lumns by a factor of 1.5-3 



Baumgartner fc Mushotzky 



20061 ). This agrees with our result 
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and indicates significant X-ray absorption due to molecular clouds rather than neutral 
hydrogen atoms, hence, the Hi column may not be a good tracer for the X-ray absorption. 

Ay, on the other hand, is caused by grains of the same heav y elements that give rise 



to X-ray absorption , therefore, it highly correlates with A^h (e.g.. 



Predehl fc Schmitt 



1995 



Giiver fc Ozel 



20091 ). Given a pulsar's p osition and distance. Ay can be estimated from the 



Drimmel et al. 



20031 ) ■ and then A^h can be deduced 



3D extinction maps of the Galaxy (e.g 
from the empirical relation A^h( cm"^)=2.21 x 10^^ Ay (mag) ( iGiiver &: Ozel 



20091). As 



shown in Figure HJ this method seems to give the best agreement between measured and 
predicted values, especially for the highest- A^h pulsars. It is worth noting that in some cases 
DMs were used to infer the pulsar distances, which then give Ay and A^h- This generally 
provides better results than directly employing the DM-A^h correlation. We believe that this 
is because the Ay map reflects the distribution of heavy elements in the Galaxy, whereas 
this crucial information cannot be obtained from DM. 



5. CONCLUSION AND OUTLOOK 

We have compiled a list of 68 pulsar A^h measurements reported in the literature 
using Chandra and XMM-Newton observations, and compared the A^h values with 
the DMs and distances. Our results show a good correlation between DM and 

with a correlation coefficient of 0.72. We obtained an empirical linear relation 
Ah (10^°cm~^) = 0.30]'^o(!,9 DM (pccm^^), implying an average ionization of 101^3%. This 
confirms the ratio of one free electron to ten neutral hydrogen atoms commonly used in 
previous studies. Our finding provides a useful tool to estimate Ah from DM. We compare 
to other Ah estimates based on the neutral hydrogen column density and Ay, and find that 
the latter gives the best results, while Hi and our empirical DM- Ah relation tend to give 
underestimates in the high- Ah regime. 



The next generation of X-ray missions, including eROSITA (iPredehl et al 



and t he proposed Neutron Star Interior Composition Explorer (NICER; 



2010) 



Gendreau et al. 



20121 ). will significantly expand the pulsar A^h sample. In addition, the foreseen Square 



Kilometer Array ( SKA) can provide parallax measurements of a few thousand radio pulsars 



( ISmits et al 



20111). Together these will allow a detailed study of the DM- Ah relation in 
different parts of the Galaxy and its connection with the Galactic structure. In addition to 
pulsars, it should be possible to compile a database of Ah measurements for other Galactic 
X-ray sources, such as stars, supernova remnants, cataclysmic variables, stellar clusters, 
white dwarfs, and X-ray binaries, and compare with their distances to build a 3D Ah map 



of our Galaxy. 
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Table 1. DM, A^h, and Distances for the 68 Pulsars Used in This Study 



PSR 


DM 




Distance* 


/ 


b 




Model° 


Ref. 




(pccm~3) 


(1020 cm-2) 


(kpc) 


(deg) 


(deg) 


(pc) 






J0030+0451 


4.333±0.001 


2.2±1.0 


n 90+0. lOp 


113.1 


-57.6 


-236 


NSAx3 


1 


J0108-1431 


2.4 ±0.2 


2 ± 2 




140.9 


-76.8 


-205 


BB 


2 


B0136+57 


73.779±0.006 


50 ±3 


9 f!+0.3p 
"^•"-0.2 


129.2 


-4.0 


-183 


PL 


3 


J0205+6449 


140.7±0.3 


41 6+°* 


3.2° 


130.7 


+3.1 


+172 


PL+RS (SNR) 


4,5 


J02 18+4232 


61.252±0.005 


8±4 


2.67<^ 


139.5 


-17.5 


-804 


PL 


6 


B0355+54 


57.1420±0.0003 


60 ±30 


l.0l°;?p 


148.2 


+0.8 


+14 


PL (PWN) 


7 


J0437-4715 


2.64476±0.00007 


^•^''-0.24 


0.156 ±0.001P 


253.4 


-42.0 


-104 


PL+NSAx2 


8 


B0531+21 (Crab) 


56.791±0.001 


32±2 


2.00° 


184.6 


-5.8 


-202 


PL 


9, 10 


J0538+2817 


39.570±0.001 


25 ±2 


1.3±0.2P 


179.7 


-1.7 


-38 


BB 


11 


B0540-69 


146.6 ± 0.2 


67 ±5 


50° 


279.7 


-31.5 


-26137 


PL (SNR) 


12, 13 


B0628-28 


34.468±0.017 


"-3 


n Q9+0-05p 
'-'■OZ_Q 04 


237.0 


-16.8 


-92.3 


PL 


14 


B0656+14 


13.977±0.013 


4.3±0.2 


0.28 ± 0.03P 


201.1 


+8.3 


+40 


PL+BBX2 


15 


J0737-3039 


48.920±0.005 


< 1 


1 1+0. 2p 
'■'^-0.1 


245.2 


-4.5 


-86 


BBx2 


16 


B0823+26 


19.454±0.004 


< 14 


n 09 + 0. 08p 
"••J^_0.05 


197.0 


+31.7 


+168 


PL 


17 


B0833-45 (Vela) 


67.99±0.01 


1 6+°-3 
'-■"-0.2 


0.28 ± 0.02P 


263.6 


-2.8 


-14 


PL (PWN) 


18 


B0950+08 


2.958±0.003 


3.2±1.3 


0.261±0.005P 


228.9 


+43.7 


+180 


PL+BB 


19 


J1016-5857 


394.2±0.2 


50±30 


8.00'^ 


284.1 


-1.9 


-263 


PL 


20 


J1023+0038 


14.325±0.010 


< 9 


1.37±0.04P 


243.5 


+45.8 


+980 


PL+NSA 


21 


J1024-0719 


6.48520±0.00008 


9+3 




251.7 


+40.5 


+318 


BB 


22 


B 1046- 58 


129.1±0.2 


90+60 


9 Q+l-2h 


287.4 


+0.6 


+29 


PL 


23 


B1055-52 


30.1±0.5 


2.7 ±0.2 


0.72d 


286.0 


+6.7 


+84 


PL+BB X 2 


15 


J1119-6127 


707.4±1.3 


200ti° 


8.4±0.4° 


292.2 


-0.5 


-78.7 


PL+BB 


24, 25 


J1124-5916 


330±2 


31±6 


r;+3h 
0-2 


292.0 


+1.8 


+153 


PL 


26 


J1231-1411 


8.090±0.001 


< 5 


0.44^ 


295.5 


+48.4 


+329 


NSA+PL 


27 


B1259-63 


146.72±0.03 


25t| 


2.3±0.4° 


304.2 


-1.0 


-40 


PL 


28, 29 


J1357-6429 


128.5±0.7 


07+20 

•^'-13 


2.50^* 


309.9 


-2.5 


-110 


PL (PWN) 


30 


J1400-6325 


563±4 


209 ± 20 


11.27<* 


310.6 


-1.6 


-313 


PL (PWN) 


31 


J 1420-6048 


358.8±0.2 


540+350 


5.61<^ 


313.5 


+0.2 


+22 


PL (PWN) 


32 


B1451-68 


8.6±0.2 


17l« 


43+0 06p 


313.9 


-8.5 


-64 


PL+BB 


33 


J1509-5850 


140.6±0.8 


2101™ 


2.62d 


320.0 


-0.6 


-28 


PL (PWN) 


34 


B1509-58 


252.5±0.3 


115±5d 


4 4+i.3h 


320.3 


-1.2 


-89 


PL (PWN) 


35 



Table 1 — Continued 



PSR 


DM 




Distance* 


/ 


& 




Moder 


Ref. 




(pc cm~^) 


(1020 cm-2) 


(kpc) 


(deg) 


(deg) 


(pc) 






J1550-5418 


830±50 


410±10 


9.55<* 


327.2 


-0.1 


-22 


PL+BB 


36 


J1614-2230 


34.4865±0.0001 




1.27=^ 


352.6 


+20.2 


+438 


BBx2 


37 


J1617-5055 


467±5 


345±20 


6.82'i 


332.5 


-0.3 


-33 


PL 


38 


J1622-4950 


820±30 


540+160 


8.73<i 


333.8 


-0.1 


-16 


BB 


39 


B1706-44 


75.69±0.05 


50±6 


9 fi+O.Bh 


343.1 


-2.7 


-122 


PL (PWN) 


40 


J1718-3718 


371.1±1.7 


130±30 


6.6<^ 


349.8 


+0.2 


+25 


BB 


41 


J1718-3825 


247.4±0.3 


70+50 
'^-13 


3.6'' 


349.0 


-0.4 


-27 


PL (PWN) 


42 


J1734-3333 


578±9 


70+40 


6.46'' 


354.8 


-0.4 


-49 


BB 


43 


J1740-5340 


71.8±0.2 


22 ±4 


2.7±0.2° 


338.2 


-12.0 


-560 


PL 


44, 45 


J1740+1000 


23.85±0.05 


10 ±2 


1.24=' 


34.0 


+20.3 


+430 


PL+BB 


46 


J1741-2054 


4.7±0.1 


15±5 


0.38<* 


6.4 


+4.9 


+33 


PL (PWN) 


47 


J1747-2809 


1133±3 


2300±150 


13.31<* 


0.9 


+0.1 


+18 


PL (PWN) 


48 


J1747-2958 


101.5±1.6 


270±10 


5° 


359.3 


-0.8 


-73 


PL (PWN) 


49, 50 


B1757-24 


289±10 


350li- 


5.22^ 


5.3 


-0.9 


-80 


PL 


51 


B1800-21 


233.99±0.05 




3.88^ 


8.4 


+0.2 


+10 


PL (PWN) 


52 


J1809-1917 


197.1±0.4 


71+60 
'^-40 


3.55=' 


11.1 


+0.1 


+5 


PL 


53 


J1809-1943 


178±5 


72±3 


3.6 ± 0.5^ 


10.7 


-0.2 


-10 


BBx3 


54 


J1819-1458 


196.0±0.4 


60±30 


3.55<* 


16.0 


+0.1 


+5 


BB 


55 


B1821-24 


120.502±0.002 


26 ±2 


5.5±0.3° 


7.8 


-5.6 


-535 


PL 


56, 57 


B1823-13 


231itl 


120tt° 


3.93'' 


18.0 


-0.7 


-47 


PL (PWN) 


58 


J1833-1034 


169.5±0.1 




4.5 ± 0.5'' 


21.5 


-0.9 


-70 


PL 


59 


B1853+01 


96.74±0.12 


120±16 


3° 


34.6 


-0.5 


-26 


VNEIX2 (SNR) 


60, 61 


J1930+1852 


308±4 


195±4 


7+3 h 
'-2 


54.1 


+0.3 


+32 


PL (PWN) 


62 


B1929+10 


3.180±0.004 


, 7+2.3 
^•'-1.7 


oi+0.09p 


47.4 


-3.9 


-21 


PL+BB 


63 


B1937+21 


71.0398±0.0002 


97 ±24 




57.5 


-0.3 


-25 


PL 


64 


B1951+32 


45.006±0.0190 


30 ±2 


3 ±2'' 


68.8 


+2.8 


+148 


PL 


65 


B1957+20 


29.1168±0.0007 


16 it 10 


2.49'' 


59.2 


-4.7 


-204 


PL 


66 


J2021+3651 


368±1 


67t? 


12.19'' 


75.2 


+0.1 


+24 


PL 


67 


J2022+3842 


429.1±0.5 


160±30 


10° 


76.9 


+1.0 


+168 


PL 


68 


J2032+4127 


114.8±0.1 


48t\l 


3.65'' 


80.2 


+1.0 


+66 


PL 


69 


J2043+2740 


21.0±0.1 


< 50 


1.8'' 


70.6 


-9.2 


-286 


BB 


17 
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Table 1 — Continued 



PSR 


DM 




Distance'' I 


b 




ModeF 


Ref. 




(pc cm"'') 




(kpc) (deg) 


(deg) 


(pc) 







J2124-3358 


4.601±0.003 


3±2 


Q or,+0.07p 


10.9 


-45.4 


-214 


PL+BB 


22 


B2224+65 


36.079±0.009 


25t« 


1.86<i 


108.6 


+6.8 


+222 


PL 


70 


J2229+6114 


204.97±0.02 




3+5o 


106.6 


+2.9 


+154 


PL 


71, 72 


J2241-5236 


11.41085±0.00003 


< 25 


0.51<^ 


337.5 


-54.9 


-417 


PL 


69 


J2302+4442 


13.762±0.006 


9+31 


1.18<i 


103.4 


-14.0 


-286 


NSA 


73 


B2334+61 


58.410±0.015 


26tf 


3.15<^ 


114.3 


+0.2 


+13 


BB 


74 



'^Distance estima tes are either parallax or Hi absorption mea surements adopted fro m 



Deller et al. 



Verbiest et al. 



1 2OI2I) and 



1 20121 ) or from the DM using the NE2001 model l|Cordes fc La2idl2002h . They are denoted by the 
letters p, h, and d, respectively. We did not attempt to quantify the uncertainties in the DM distances. The 
designation "o" indicates distance estimates based on other arguments; see the references for details. 

''Vertical distance from the Galactic Plane calculated using the source distance and Galactic latitude b. 

'^Spectral models used to obtain Nn: blackbody (BB), power law (PL), and neutron star atmosphere (NSA). A^h 
values determined from the associated SNRs or PWNe are noted. The SNR spectra were fitted with Raymond- 
Smith (RS) and non-equilibrium ionization (VNEI) models. 

I 2OI0I) reported A^h — 1-15 X 10 cm , but mentioned that it is consistent with the previous 
(2002), which gave A^h = 9.5 X 10^^ cm~^. Therefore, the former is assumed to be a 



Schock et al 



result from 



Gaensler et al 



typographical error and we adopt the value 1.15 X 10^^ cm ^ . 

Note. — Uncertainties and upper limits on Ajj are all scaled to the 90% confidence level, i.e. 1.6(7. 
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Fig. 1. — A^H versus DM for 68 pulsars, (a) Data points in gray color, including the Vela 
pulsar, PSR B0540— 69, and A^h measurements with fractional uncertainties larger than 80% 
or only upper limits are not used in the fit. The red solid and dotted lines show the best 
linear fit with the 90% confidence interval, and the blue dashed line is the same fit ignoring 
measurement uncertainties. These correspond to 10^3% and 4% ionization, respectively. 
Uncertainties in DM are negligible. The same plot is shown in (b) and (c) with different 
color schemes, indicating the pulsar vertical distance from the Galactic Plane and their 
Galactic longitudes, respectively. 
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Fig. 2. — A^H versus distance for Galactic pulsars. PSR B0540— 69 is not shown here. 

(a) The colors indicate different types of distance measurement. We did not estimate the 
uncertainties for the DM distances. The same plot with different color schemes is shown in 

(b) and (c), indicating the pulsar vertical distance from the Galactic Plane and their Galactic 
longitudes, respectively. 
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Fig. 3. — Same as Figure Ej but for DM versus distance. Uncertainties in DM are not plotted, 
as they are smaller than the data points. 
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Fig. 4. — Comparison between iVn estimates and measurements. The underlying plots in 
(a) and (b) are the same as Figures [1] and [21 respectively. The red crosses indicate the total 
actic Hi column densities for each pulsar, given by 21cm radio observations 
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20051 ). The blue triangles show estimates based on Ay (see text). In (c), the 



open circles represent predictions from our best-fit DM-A^h relation. 



